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ABSTRACT 

Wc present Hubble Space Telescope Observations of (596) Scheila during its 
recent dust outburst. The nucleus remained point-like with absolute magnitude 
V{1, 1, 0) = 8.85±0.02 in our data, equal to the pre-outburst value, with no sec- 
ondary fragments of diameter >100 m (for assumed albedos 0.04). We find a 
coma having a peak scattering cross-section ~2.2xl0^ km^, corresponding to a 
mass in micron-sized particles of ~4xl0^ kg. The particles are deflected by so- 
lar radiation pressure on projected spatial scales ~2xl0^ km, in the sunward 
direction, and swept from the vicinity of the nucleus on timescales of weeks. The 
coma fades by ~30% between observations on UT 2010 December 27 and 2011 
January 04. The observed mass loss is inconsistent with an origin either by rota- 
tional instability of the nucleus or by electrostatic ejection of regolith charged by 
sunlight. Dust ejection could be caused by the sudden but unexplained exposure 
of buried ice. However, the data are most simply explained by the impact, at 
~5 km s~^, of a previously unknown asteroid ~35 m in diameter. 

Subject headings: minor planets, asteroids: individual ((596) Scheila)); comets; 
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1. Introduction 

Main belt asteroid (596) Sclieila (formerly 1906 UA and hereafter "Scheila") was dis- 
covered photographically in 1906 by August Kopff. It is a large body, with equivalent 
circular diameter D = 113±2 km and a visual geometric albedo = 0.038±0.004 (Tedesco 
et al. 2002). The normalized (to 5500A) optical reflection spectrum of Scheila (Bus and 
Binzel 2004) has a gradient S' = 6.2±1%/1000A (the uncertainty is our estimate of the sys- 
tematic error; the statistical error is only ±0.04%/1000A). The low albedo and the spectral 
slope together suggest that Scheila can be spectrally classified as a primitive P- or D- type 
(Dahlgren and Lagerkvist 1995). Under the assumption of a bulk density p = 2000 kg m~^ 
and a spherical shape, the approximate escape speed from Scheila is Ve = 60 m s~^. 

Scheila's orbit lies in the outer asteroid belt (semimajor axis, a = 2.926 AU, eccentricity, 
e = 0.1644, and inclination, i = 14.7°). The Tisserand parameter with respect to Jupiter, 
Tj = 3.21, is typical of asteroids and lies far above the Tj = 3 dynamical dividing line 
separating comets from asteroids (Kresak 1980). However, Larson (2010) discovered a comet- 
like appearance in observations taken UT 2010 December 11.44-11.47 with the 0.68 meter 
Catalina Schmidt telescope. Observations with the same telescope on UT 2010 Dec 3.4 
showed a slightly diffuse appearance and an integrated magnitude V = 13.2, about 1.3 mag 
brighter than in the previous month, when the object appeared point-like (Alex Gibbs, 
reported in Larson 2010). With a main-belt orbit (Tj > 3) and a comet-like morphology, 
Scheila satisfies the definition of a main-belt comet (MBC: Hsieh and Jewitt 2006). It is the 
seventh known example (and by far the largest) of this class (Figure [T]). The MBCs appear 
to be of diverse origins, including examples likely to be driven by the sublimation of near 
surface ice (e.g., 133P/Elst-Pizarro; Hsieh and Jewitt 2006) and others likely to result from 
recent, probably colhsional, disruption (e.g., P/2010 A2; Jewitt et al. 2010, Snodgrass et al. 
2010). 

Here we report initial observations with the Hubble Space Telescope (HST), taken to 
examine Scheila at high angular resolution shortly after its outburst. 

2. Observations 

We secured two orbits of HST Director's Discretionary Time on UT 2010 Dec. 27.9 and 
UT 2011 Jan 04.9 with the WFC3 camera (Dressel et al. 2010). On each orbit we took six 
exposures each of 4 s duration to image the near-nucleus environment and four exposures 
each of 390 s duration to examine the low surface brightness coma. The short integrations 
employed the F621M filter (central wavelength ~ 6210A and FWHM [full width at half 
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maximum] ~ 640 A) to prevent saturation while maintaining the point-spread function of 
the image. The long integrations used the F606W filter (Ac ~ 6OOOA and FWHM ~ 2300 A) 
to maximize sensitivity to low surface brightness coma. At the time of observation the 0.04" 
pixels corresponded to 66 km at the comet, so that the Nyquist sampled (2 pixel) resolution 
was 133 km. A brief journal of observations is given in Table [TJ 

The images from each orbit were combined into composites having 24 s effective inte- 
gration for the nucleus and 1560 s integration for the coma. The nucleus (A^ in Figure |2]) 
was found to be point-like, with FWHM = 2.2 pixels (0.09"). Photometry of the nucleus 
was obtained within a circular projected aperture of radius 0.4", with sky subtraction from 
a concentric annulus with inner and outer radii 0.6" to 0.8", respectively (Table [2]). The 
absolute magnitude, \^(1,1,0), (the magnitude if observed from heliocentric and geocentric 
distances R = A = 1 AU and at phase angle a = 0°) is computed from 

y(l,l,0) = y-51ogio(i?A)-/(a). (1) 

in which / (a) is the phase function. For the latter we adopt the H — G formalism of Bowell et 
al. (1989), with G = 0.076±0.060 as found from pre-outburst observations (Warner 2010). 
As seen in Table |2| the nucleus apparent magnitude brightened by ~0.12 mag in the 8 
day interval between the HST observations, but this brightening is consistent with the 
changing observing geometry of the object, since V^(1,1,0) = 8.85±0.02 remains constant. 
The derived absolute magnitudes from HST are consistent with the pre-outburst value, 
8.84±0.04 (Warner 2010), showing that near-nucleus dust is insignificant. 

On larger scales, the coma appeared diffuse and asymmetric, being brighter and more 
extended to the North than to the South of the nucleus on both visits {A and B in Figure 
|2|. The sunward extension of the coma was ~10", corresponding to ~ 16, 000 km in the plane 
of the sky. The overall form of the coma suggests that dust particles launched sunwards 
are being slowed by solar radiation pressure and so are concentrated near their turn-around 
points, giving rise to a bright-edged parabola (Figure |2]). We connect s, the turn-around 
distance along the Sun-comet line to u, the initial sunward particle speed from 

= 2/3gQS (2) 

where (3 is the dimensionless radiation pressure factor and Qq is the gravitational acceleration 
towards the Sun. (Strictly, the data provide only a lower limit to s because of the effects of 
projection (the phase angle was ~13°, c.f. Table [T|). Substituting for gQ we obtain 
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where G is the gravitational constant, R the hehocentric distance and Mq the mass of the 
Sun. The characteristic ejection velocity must be m > Ve, since slower particles should fall 
back to the nucleus under gravity. Substituting -u > 60 m s~^ in Equation [s] we obtain /3 > 
0.2. The magnitude of /3 is inversely related to particle size; (3 > 0.2 is compatible with 
dielectric particles having radii ~0.1 to 1 /im (Bohren and Huffmann 1983). We conclude 
that the diffuse coma of Scheila is populated by small dust grains just like those in other 
comets at 3 AU. In this model, the characteristic travel time from the nucleus to the apex 
of motion is r = 2s/u, or r = 5x10^ s (~1 week), while the residence time in the HST field 
of view would be a month or more. We note that the rapid fading of the coma on timescales 
of approximately a month strongly suggests that we are not observing large (>10 fim) slow- 
moving grains in Scheila, in stark contrast to the situation for P/2010 A2 (Jewitt et al. 
2010). 

To measure the brightness of the surrounding coma, we first digitally removed field star 
and galaxy image trails that were not already cancelled by the image combination process. 
For this purpose we replaced afflicted pixels with the average pixel value measured in a 
surrounding region. The coma brightness was measured using a projected circular aperture 
64" in radius, with sky subtraction from the median signal measured within a contiguous 
aperture having outer radius 80". The integrated magnitudes are brighter than the nucleus 
by about AV = 1.3 mag on Dec 27 but only by AV = 1.0 mag on Jan 04 (cf. Table [2]). 
This ~30% fading of the coma in the 8 days between measurements is broadly consistent 
with the timescale deduced above from radiation pressure sweeping, assuming the coma is 
not being continuously replenished. 

We calculated the effective scattering cross-section of the coma from Cc = C„(10°'^^^ ~ 1); 
where C„ = vrr,^ = 1.0x10"^ km^ is the geometric cross-section of the nucleus. We further 
calculated the mass of dust particles in the coma from the scattering cross-section using 



M, = paC, (4) 

where p is the particle density, taken to be p = 2000 kg m~^, and a is the average particle 
radius in the coma. We take a = 1 /im, corresponding to the upper limit estimated above from 
radiation pressure effects. This size is consistent with the radiation pressure considerations 
discussed above and is also characteristic of optical observations of normal comets since much 
smaller particles are inefflcient scatterers of optical photons while much larger particles are 
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rare. The resulting dust cross-sections and masses are listed in Table |2} The computed mass 
in micron-sized grains is strictly a lower limit to the total mass, since large particles may 
hold significant mass while presenting negligible scattering cross-section. Both Cc and Mc 
decreased between December 27 and January 4 by about 30%, indicating that the escape 
of particles from the projected 64" radius photometry aperture substantially exceeded the 
supply of fresh particles from the nucleus in this period. 

In addition to the diffuse coma, an approximately linear tail (the "spike") is evident 
(marked C in Figure |2]), with a position angle 277° ± 1° on Dec 27 and 273° ± 1° on Jan 04. 
The position angles are slightly different from both the antisolar direction (position angles 
276° and 269° on Dec 27 and Jan 04, respectively) and the projected orbit (position angle 
286° on both dates). Therefore, it is not possible to interpret the spike as a simple synchrone, 
as was done for P/2010 A2. Instead, it must consist of particles whose motion is determined 
by their initial velocity as well as radiation pressure. 

We set upper limits to the allowable brightness of co-moving companions to Scheila by 
digitally adding scaled versions of the (unsaturated) nucleus. At projected distances >6" 
(10,000 km), any companion with \^ < 28 could not escape detection. Assuming that any 
companion has the same albedo as Scheila, a limit to the diameter of any secondary is set 
at < 0.1 km. 



3. Discussion 

Two mechanisms considered previously as possible explanations for activity in main-belt 
comets are inoperable on Scheila, as a result of its large size and slow rotation. First, mass loss 
through rotational instability is ruled out by the measured rotational period, P = 15.848 hr 
(Warner 2006). The latter greatly exceeds the critical period at which the centripetal ac- 
celeration at the surface of a sphere equals the gravitational acceleration, assuming a bulk 
density p = 2000 kg m~^. Second, the ejection of grains through electrostatic charging of the 
surface can be ruled out since the speeds generated electrostatically (v ~1 m s~^) (Rennil- 
son and Criswell 1974) are far smaller than the ~60 m s~^ escape speed from the nucleus. 
Another mechanism must be responsible for the ejection of dust. 

The simplest explanation is that Scheila ejected material after being struck by another, 
much smaller, asteroid. This is the explanation proposed elsewhere for the inner-belt MBC 
P/2010 A2 (Jewitt et al. 2010). In this interpretation, the estimated dust mass (Table |2]) 
gives a crude estimate of the impactor properties. The velocity dispersion amongst main- 
belt asteroids is Vi ~5 km s~^, about 10^ times the escape velocity from Scheila. At such a 
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high speed, experiments show that a projectile can excavate 10^ times its own mass from 
the target. However, the bulk of the ejecta moves too slowly to escape. Specifically, the 
mass of the ejecta leaving with Ve/vi >10~^ is comparable to the impactor mass (Figure 4 
of Housen and Holsapple 2011). Therefore, the impact hypothesis requires that Scheila be 
struck by a body having mass M ~4.4xl0'' kg (Table [2]). If also of density p = 2000 kg m~^, 
the impactor diameter would have been (ij ~ 35 m and the kinetic energy of the impact 
~5.5xl0^^ J (about 0.1 MTonnes TNT equivalent). The resulting crater on Scheila would 
be of diameter ~ 400 m. 

Another possibility is that the coma is produced by sublimation of surface ice, as in 
133P/Elst-Pizarro (Hsieh and Jewitt 2006). Radio spectral line observations limit the pro- 
duction rate from Scheila to logioQo/f < 26.3 in the period December 14 - January 04, 
corresponding to mass production rates <6 kg s~^ in water (Howell and Lovell 2011). How- 
ever, these production limits are model-dependent, and cannot exclude the possibility that 
Scheila outgassed more strongly at earlier times, launching dust from the surface by the 
action of gas drag forces. 

The energy balance equation for sublimating ice may be written 
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in which Fq = 1360 W m~^ is the Solar Constant, Rau is the heliocentric distance expressed 
in AU, r„ (m) is the nucleus radius, A is the Bond albedo, e is the emissivity, T (K) 
the effective temperature of the nucleus, L(T) (J kg^^) the temperature-dependent latent 
heat of sublimation and dm/dt (kg m~^ s~^) the rate of sublimation of the ice per unit area. 
Parameter x iii Equation ([s]) is a proxy for the surface temperature variation over the surface 
of the nucleus, itself dependent on the nucleus shape, thermal properties, rotation period and 
spin axis direction relative to the Sun. To solve Equation ([s]) requires additional knowledge 
of L{T) and dm/dt. For this, we use the saturation vapor pressures and latent heat data for 
ice from Washburn (1926). 

Allowable values of x li^ in the range 1 < x ^ 4. The high temperature limit, x = 1; 
describes a flat plate oriented with its normal pointing towards the Sun, because then the 
absorbing and radiating areas are identical. This would approximate, for instance, the 
heating of the sunward pole on a nucleus whose spin-vector points towards the Sun. The 
cold limit, x = 4, corresponds to an isothermal, spherical nucleus in which solar power is 
absorbed over an area vrr^ but radiated from 47rr^. True cometary nuclei probably possess 
intermediate effective values of x- For example, a spherical nucleus with the Sun on the 
polar axis corresponds to x = 2, while a nucleus with the Sun in the equator and with 
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rotation so rapid as to maintain temperature along lines of constant latitude has x = ^) 
and so on. The point is that the rotational and thermophysical properties of a nucleus are 
usually unknown, so that the effective value of Xi ^"^^ thus the sublimation rate, cannot be 
accurately calculated. 

We used Equation (|5| to calculate dm/dt as a function of Rau for spherical model comets 
with water ice exposed at the surface. We used A = 0.03 and assumed e = 0.9. The limiting 
values X = 1 ^iiid x = 4 were assumed in order to bracket the range of likely behaviors to 
be expected from the real comets. Results are shown in Figure ^ for 2 < Rau < 4 AU. 
Also shown in Figure ^ is the critical radius, Oc, defined as the radius of the largest dust 
particle that can be ejected from Scheila against gravitational attraction to the nucleus, 
given sublimation at the rate computed from Equation 

Figure (|3| shows that allowable temperature regimes on Scheila permit more than two or- 
ders of magnitude variation in the sublimation rate per unit area, from ~3xl0~^ kg m~^ s~^ 
in the hot case to ~8xl0~^ kg m~^ s~^ in the cold case, at Rau = 3 AU. If we assume that 
the dust/gas ratio of ejected material is of order unity, and that the dust ejection occurred 
over a period of less than 1 month (2.5x10^ s), then the area of exposed ice needed to sup- 
ply the M ~2xl0^ kg coma mass is A = 0.3 km^ in the high temperature case rising to 
A = 100 km^ in the low temperature case. The surface area of Scheila is ~40,000 km^, so 
that even sublimation from the larger area corresponds to only 0.25% of the surface. 

The critical dust radius for ejection at 3 AU varies from about 10~^ m in the hot case 
to 6xl0~^ m in the cold case. Thus, at this distance from the Sun, gas sublimated from 
exposed water ice is capable of launching dust particles above the escape speed no matter 
what the temperature regime on the surface. Even at aphelion and with the lowest plausible 
surface temperatures, sublimated ice can eject grains with ac ~ 0.1 /im, the minimum size 
for efficient scattering of optical photons. Therefore, both in terms of the area of exposed 
ice required, and in terms of the particle speeds induced by gas drag, ice sublimation is a 
plausible driver of the activity observed in Scheila at 3 AU, and potentially could supply 
observable coma at any position in the orbit. 

The ice hypothesis has two problems, however. First, surface ice on Scheila is long-term 
unstable to sublimation and must be protected by a refractory mantle (Schorghofer 2008). 
What process exposed the ice in early December? The most obvious possible trigger is an 
impact which, by the discussion above, would need to expose an area of ice A > 0.3 km^ 
(corresponding to a circle 300 m in radius or larger). But, as discussed above, an impact this 
large could produce an impulsive coma as observed even without the need for a sublimating 
ice source. 
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Second, if ice sublimation is the driver of the activity, why does the coma fade on 
timescalcs of r ~ 8 days (7x10^ s)? At dm/dt = SxlO^'"* kg m^^ s^^, ice would recede into 
the surface by a distance S = Tp~^{dm/dt) which, by substitution gives S ~ 0.01 m (1 cm). 
This is smaller than the diurnal thermal skin depth (d = 0.07 m [7 cm]) making it hard to 
see why the coma formation, once started, would so soon stop. For these two reasons, we 
conclude that the impact hypothesis provides the most simple and most likely explanation 
of the observed activity. Scheila now stands, after P/2010 A2, as the second example of an 
impact-activated main-belt comet. 

4. Summary 

From images of (596) Scheila obtained using the Hubble Space Telescope we find that: 

1. The coma has effective peak cross-section 2.2x lO'' km^, is shaped by radiation pressure, 
and fades on a timescale consistent with radiation pressure clearing. 

2. The coma dust particles have radiation pressure factors /3 >0.2, corresponding to 
micron-sized dielectric particles. The mass in micron-sized dust particles is 4x10^ kg. 

3. The scattering properties of the nucleus appear unchanged by the ejection of dust. 

4. No near-nucleus fragments or structures are apparent in the Hubble data. 

5. The measurements are consistent with dust ejection by impact into Scheila of a ~35 m 
diameter projectile. 

Based on observations made with the NASA/ESA Hubble Space Telescope, with data 
obtained from the archive at the Space Telescope Science Institute (STScI). STScI is operated 
by the association of Universities for Research in Astronomy, Inc. under NASA contract 
NAS 5-26555. We thank the Director's Office at the STScI for granting us Discretionary Time 
to make these observations and Alison Vick, Larry Petro, and other members of the STScI 
ground system team for their expert help in planning and scheduling these observations. Bin 
Yang and an anonymous referee supplied useful comments on the manuscript. DJ thanks 
Mike Hicks of JPL for supplying his Palomar observations. Support is appreciated from 
NASA's Planetary Astronomy program to DJ. 
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Table 1. Journal of Observations 



Instrument 












WFC3/F606W/F621N 


2010 Dec 27.9 


3.085 


2.338 


13.7 


-2.6 


WFC3/F606W/F621N 


2011 Jan 04.9 


3.073 


2.254 


11.9 


-3.5 



*UT Date of the observation 
''Heliocentric distance in AU 
'^Geocentric distance in AU 
^Phase (sun-object-Earth) angle in degrees 
^Out-of-plane angle in degrees 
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Table 2. Photometry 



Date 


Vn" 


^(1,1,0)" 










2010 Dec 27.9 

2011 Jan 04.9 


13.98±0.02 
13.86±0.02 


8.84±0.02 
8.86±0.02 


12.63 
12.86 


1.26 
1.00 


2.2x10^ 
1.5x10^ 


4.4x10^ 
3.0x10^ 



^Nucleus y-band magnitude measured within a 0.4" radius aperture. The 
apparent V-band magnitude was computed from the observed count rate "C" 
in electrons s~^) using V = -2.5 log C + Z, where Z = 24.45 for the F621M 
fiher and 25.99 for the F606W filter (Kalirai et al. 2009). 

^Absolute V magnitude of the nucleus (Equation [l| 

■^Total y-band magnitude measured within a 64" radius aperture 

'^AV = Vn-V 

"^Effective scattering cross-section in km^ 
^Effective dust mass, in kg 
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Semimajor Axis [AU] 



Fig. 1. — Distribution of the known MBCs (red circles) in the semimajor axis vs. orbital 
eccentricity plane. The corresponding distributions of asteroids (orange circles) and comets 
(blue circles) are shown for comparison. Objects above the diagonal arcs cross either the 
aphelion distance of Mars or the perihelion distance of Jupiter, as marked. The semimajor 
axes of the orbits of Mars and Jupiter are shown for reference, as is the location of the 2:1 
mean-motion resonance with Jupiter. 



Fig. 2. — Composite images from UT 2010 Dec 27 (top) and UT 2011 Jan 4 (bottom), each 
of 1560 s exposure. The images have North to the top, East to the left and are shown with 
identical stretches. In the top panel, N marks the nucleus, A and B are the north and south 
extensions of the small particle coma, and C is the spike discussed in the text. Marked 
arrows show (AS) the projected anti-solar direction and (O) the projected orbit. The gap 
between chips in the detector is marked. 
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Heliocentric Distance, R [AU] 

Fig. 3. — Two-axis plot showing, on the left, the equilibrium water ice sublimation rate 
(solid lines) and, on the right, the maximum ejectable dust grain radius (dashed lines), 
both as functions of the heliocentric distance. Red(black) curves correspond to the mini- 
mum(maximum) temperature extremes, as discussed in the text. The shaded region marks 
the range of heliocentric distances swept by Scheila each orbit. 



